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Investigations of the optical response of subwavelength structure arrays milled into thin metal 
films has revealed surprising phenomena including reports of unexpectedly high transmission of light. 
Many studies have interpreted the optical coupling to the surface in terms of the resonant excitation 
of surface plasmon polaritons (SPPs), but other approaches involving composite diffraction of surface 
evanescent waves (CDEW) have also been proposed. We present here a series of measurements on 
very simple one-dimensional (1-D) subwavelength structures with the aim of testing key properties 
of the surface waves and comparing them to the CDEW and SPP models. 



I. INTRODUCTION 

Initial reports of dramatically enhanced transmission through arrays of subwavelength holes in thin films and 
membranes [1-3] have focused attention on the physics underlying this surprising optical response. Since the early 
experiments were carried out on metal films, surface plasmon polaritons [4, 5] were invoked to explain the anomalously 
high transmission and to suggest new types of photonic devices [5] . Other interpretations based on "dynamical 
diffraction" in periodic slit and hole arrays [6, 7] or various kinds of resonant cavity modes in 1-D slits and slit arrays [8, 
9] have also been proposed. Reassessment of the earlier data by new numerical studies [10] and new measurements [11] 
have prompted a sharp downward revision of the enhanced transmission factor from ~ 1000 to ~ 10 and have 
motivated the development of a new model of surface wave excitation termed the composite diffracted evanescent 
wave (CDEW) model [11]. This model builds a composite surface wave from the large distribution of diffracted 
evanescent modes (the inhomogeneous modes of the "angular spectrum representation" of wave fields [12]) generated 
by a subwavelength feature such as a hole, slit, or groove when subjected to an external source of propagating wave 
excitation. The CDEW model predicts three specific surface wave properties. First, the surface wave is a composite 
or "wave packet" of modes each evanescent in the direction normal to the surface. The surface wave packet exhibits 
well-defined nodal positions spaced by a characteristic wavelength, Asurf ; second, the appearance of the first node at 
a distance of Asurf /2 from the subwavelength launch site (essentially a phase delay of 7r/2 with respect to the E-field 
of the external driving source); and third, an amplitude decreasing inversely with distance from the launch site. We 
present here the results of a series of experiments on very simple 1-D subwavelength surface structures designed to 
investigate these predictions and thus assess the validity of the model. 

II. SUMMARY OF THE CDEW MODEL 

The essential elements of the CDEW model can best be summarised with reference to Fig. 1. It is based on a 
solution to the 2-D Helmholtz equation in the near field and subject to the slab-like boundary conditions of a slit in 
an opaque screen. The basic expression describing the scalar wave is 

[V^ + e]E{x,z)=0 (1) 

with = d^/dx^ -f d^/dz^, k — 27r/A and E{x,z) the amplitude of the wave propagating in the x,z directions. 
Kowarz [13] has written down the solution to this equation for the case of an incident plane wave propagating in air 
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FIG. 1: Essential elements of the CDEW model. The incoming plane wave Ei with ko — 2n/Xo in air {n = 1) is linearly 
polarised parallel to the plane of the structure and perpendicular to the slit of subwavelength width d. A fraction of the 
incoming light -Esurf forms the composite diffracted wave in the ±x directions, and the blue trace (displaced above the surface 
for clarity) shows Esurf- The alternating blue and red loops indicate the field lines induced by the CDEW near the surface of 
the silver film. The green trace (offset above the blue trace for clarity) shows the cosine representation of the CDEW expressed 
by Eq. 4 and closely approximating Eq. 2a for > 3/4 A. The wavelength Asurf = Ao/riaurf where risurf is the surface index of 
refraction. 



(n = 1) with amplitude Ei and propagation vector fco impinging on a slit of width d in an opaque screen. Specifying 
the coordinates as shown in Fig. 1, the field solution E^v for the modes evanescent in z at the z = boundary is 



E- 

E^^{x,z = 0) = ^<jSi 



k[ x+- 



Si 



for |a;| > d/2 



E, 



TT — Si 



k{x + 



Si 



k X 



for \x\ < d/2 with Si(a) = / 

Jo 



" sin(t) 



(2a) 



dt (2b) 



The kz evanescent modes are determined by a conservation-of-energy criterion, 
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fr2 



kx > kn 



(3) 



The form of the inhomogeneous or evanescent field on the z = boundary is shown in Fig. 1. At transverse displace- 
ments from the slit \x\ > d/2, the evanescent component of the field at the surface Ecv{x, z = 0) can be represented 
to good approximation by the expression 



Ei d . . , 

Ecv — cos (fcsurfa; + 7r/2) 

TT X 



(4) 



that describes a damped wave with amplitude decreasing as the inverse of the distance from the launching edge of 
the slit, a phase shift 7r/2 with respect to the propagating plane wave at the midpoint of the slit and a wave vector 
fcsurf = 27r/ Asurf- The wavelength of the CDEW on the surface Asurf = Aq /risurf where Uguri is the surface index of 
refraction (empirically, risurf — 1.04). This surface wave is actually a composite superposition of k^ modes evanescent 
in z, with \kx\ > fcg and directed along the ±x axes. 
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Equation 5 generalises the expressions of Eqs. 2a, 2b to include the evanescent components above the z = plane. 
When the composite evanescent wave encounters a surface discontinuity (a slit for example), a fraction of the surface 
wave is reconverted to a distribution of "homogeneous" or propagating modes |k| = 2tt/Xq at the site of the slit. In a 
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FIG. 2: The incoming plane wave Ei impinges on the subwavelength sht (or hole) and a groove milled on the input side. The 
evanescent -Esurf wave originates on the surface at a slit-groove distance Xsg and is indicated in blue. In the model proposed 
in [11] CDEWs travel along the surface toward the slit where they reconvert to a propagating field E^i and interfere with Et, 
the propagating field directly transmitted through the slit or hole. The superposed output field Eo = Et + Eai propagates into 
the 2 > half-space and the intensity of the interference figure I{9) is detected in the far field. 



laser source 
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FIG. 3: Goniometer setup for measuring far-field light intensity and angular distributions. A stabilised single mode CW diode 
laser, locked to a wavelength of 852 nm and modulated at 850 Hz by a chopper wheel, is injected into a single-mode fibre and 
focused onto the nanostructures mounted in a x-y translation stage as shown. A stepper motor drives the goniometer arm, and 
the chopped light intensity detected by the photodiode is fed to a lock-in amplifier. Output from the lock-in is registered by 
the PC that also drives the stepper motor. For the input-side experiments described here the detector was always positioned 
at e = 0° . 

practical experiment, any real planar structure has two surfaces: an "input side" in the half-space z < 0, containing 
the incoming plane wave, and an "output side" in the half-space z > 0, containing the far- field propagating modes 
issuing from the output surface and a photodetector. Experiments can be carried out by fabricating subwavelength 
grooves on the input side, the output side or both. The measurements reported here concern only the input-side 
experiments (Fig. 2). Results for output-side experiments will be reported later. 

III. MEASUREMENTS AND RESULTS 

Measurements of the optical response of the slit-groove and hole-groove structures were carried out using a home- 
built goniometer shown in Fig. 3, details of which are described in the caption of Fig. 3 and in the Methods section. 
We have carried out a series of measurements on simple 1-D structures to test the "signature" predictions of the 
CDEW model, viz. (1) a composite surface wave expressed by Eq. 5 and approximately represented by a damped 



FIG. 4: Scaiming electron microscope (SEM) image of 
one of the series of singlc-slit, single-groove structures FIB 
milled into a 400 nm thick silver layer deposited on flat 
quartz microscope slides 1 mm thick. The width of both 
the slit and the groove is 100 nm, the height 20/irn and 
the groove depth ~ 100 nm. The distance Np is the 
pitch increment p = 104 nm multiplied by the number of 
increments N. 



FIG. 5: SEM image of one of the series of single-groove, 

single-hole structures fabricated similarly to the single- 
groove, singlc-slit structures of Fig. 4. The silver layer for 
the groovc-holc structures is 260 nm thick, the width and 
depth of the groove is 100 nm and 70 mn respectively, 
and the diameter of the hole is 300 nm. The distance 
Np is the pitch increment p = 104 nm multiplied by the 
number of increments N 



wave, Eq. 4; (2) a phase shift of 7r/2 between the CDEW and the driving source plane wave and (3) a wave amplitude 
that decreases inversely with distance from the launching groove. Figures 4, 5 show one of the series of structures 
consisting of one slit and one groove and one hole and one groove, respectively. The slit-groove distance Xsg or 
hole-groove distance Xhg is indicated as Np where p is the basic unit of distance increment, the "pitch," and N is the 
number of increments. The pitch p was taken to be 104 nm, approximately one-eighth the wavelength of the surface 
wave and N was varied from 4 to 59. Structural details of these devices are described in the captions of Figs. 4, 5 
and in the Methods section. The slit(hole)-groove structures were mounted facing the input side and exposed to 
plane-wave radiation from the focused TEMqo laser source. Measurements of light intensity on the output side in the 
far field, 200 mm from the plane of the structures, were carried out on the slit-groove structures using the goniometer 
setup described in the Methods section. The results are shown in Figs. 6, 7. They show an oscillatory fringe pattern 
with amplitude damping out to a distance of ~ 3 — 4 /im and maintaining an essentially constant amplitude from that 
point out to the distance limit of the measurements. As indicated in Fig. 2, the fringe pattern results from interference 
between the mode directly propagating through the slit (hole) at the input side Et and a surface wave originating 
from the single-groove structures -Bgurf • The wave -Esurf is reconverted to a propagating mode at the slit or hole, and 
it is this propagating mode that interferes with . The frequency and phase of the interference pattern is a function 
of the slit (hole)-groove optical path and any intrinsic phase shift of the surface wave itself. The normalised intensity 
I/Iq of the superposition term is given by 



where a = Eguii /Ei is the fractional amplitude of the surface wave launched from the incoming field Ei at the groove 
site, and /? is the further fraction of this surface wave reconverted to a propagating wave in the slit, E^i = pEsuii = 
(3aEi. The fractional amplitude of the directly transmitted component Et is 5 and the phase difference 7i between 
Et and Egi is the sum of two terms, 



The first term ksuiiXsg{hg) is the phase accumulated by the surface wave propagating from the groove to the slit (hole) 
and the second term (fi^t is any phase shift intrinsic to the surface wave. The term ipint includes the "signature" 
shift of the CDEW plus any phase shift associated with the groove width and depth. Figures 6, 7 present a direct 
measure of the normalised amplitude damping with distance, rji = rji {x) and the period and phase of the oscillations, 
from which the wavelength Agurf of the surface wave, the phase ipi^t, and the effective surface index of refraction 
nsurf can be determined. Analysis of the frequency spectrum of the fringe pattern for the slit (hole) structures results 
in the determination of a surface wavelength Agurf = 819(811) ± 8nm and an effective surface index of refraction 
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FIG. 6; Normalised far-field intensity ///q as a function of 
slit-groove distance Xsg for series of single-slit, single-groove 
structures mounted facing the input side with respect to 
plane wave excitation. Points are the measured data through 
which the solid line, Eq. 8a, is fitted with parameters Hsi ~ 
0.13±0.01, Ksi = 0.12±0.01/imandv3fnt = 0.81±0.027r. Error 
bars were determined from variations in the measured inten- 
sities of the six nominally identical naked slits (no flanking 
groove) used for normalisation of each measurement. Anal- 
ysis of the frequency spectrum of the fringe pattern for the 
slit-groove structures results in the determination of a surface 
wavelength Asurf = 819 ± 8nm and an effective surface index 
of refraction rigurf = 1.04 ± 0.01. 
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FIG. 7: Normalised far-field intensity I/Iq as a function 
of hole-groove distance Xhg for series of single-hole, single- 
groove structures mounted facing the input side with re- 
spect to plane-wave excitation. Points are measured data 
through which the solid line, Eq. 8b, is fitted with parame- 
ters Hhi = 0.10 ± 0.02, KM = 0.08 ± 0.03 ^im, and ipilt = 
0.55 ± 0.05 TT. Error bars were determined from variations in 
the measured intensities of the six nominally identical naked 
holes (no flanking groove) used for normalisation of each mea- 
surement. Analysis of the frequency spectrum of the fringe 
pattern for the hole-groove structures results in the determi- 
nation of a surface wavelength Asmf = 811 ± 8nm and an 
effective surface index of refraction riaurf — 1.05 ± 0.01. 



n-aurf = 1.04(1.05) ± 0.01. The amplitude rji of the oscillatory term depends on the slit(hole)-groove distance, and 
Figs. 6, 7 show that rji falls of with increasing distance. This fall-off is fit to an expression with two terms: an inverse 
distance dependence term plus a constant term. 

r]f{xsg) cos(7i) = ( — + ^isi) cos{ksuriXsg + (fitit) (8a) 

Vi'{xhg)cOs{-fi) = ( — + fi hi) COs{ksurf Xhg + tPint) (8b) 
\ ^hg J 

The best-fit values for fJ,,K,(p are indicated in the captions of Figs. 6,7 and in Table I for slit and hole structures, 
respectively. The subscript i and superscripts si, hi on 77 refer to input-side, slit and hole measurements, respectively. 



IV. DISCUSSIONS AND CONCLUSIONS 



The measured interference fringes on slit (hole) structures exhibit the presence of a surface wave with wavelength 
819(811) ± 8 nm and therefore a surface index of refraction rigurf = 1.04(1.05) ± 0.01. The amplitude behavior of 
these fringes is also similar. Both the slit-groove and hole-groove structures exhibit an initial amplitude fall-off with 
increasing distance, damping to an essentially constant amplitude at a distance ~ 6 /im. This behavior is fit to the 
expressions in Eqs. 8a, 8b. The fitting procedure is a linear regression varying the relative contributions of the constant 
and decaying amplitude terms, ^isi(hi) ^^^id Ksi{hi) ^^id the intrinsic phase shift (/j^J^'''-*. The results are summarized in 
Table I. 

How do these results compare to CDEW or SPP models? In the CDEW picture, the groove launches a surface wave 
on the input side of the silver film that is detected by interference with the directly transmitted wave through the hole 
or slit, in the far field, on the output side of the structure. The amplitude of this surface wave is predicted to damp 
as the inverse distance between the groove and the slit or hole. Figures 6, 7 show an initial decrease in amplitude 
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TABLE I: Fit parameters and measured Ug for slit and hole structures. Error bars were determined from variations in the 
measured intensities of the six nominally identical naked slits and holes (no flanking groove) used for normalisation of each 
measurement in Figs. 6, 7. 



parameter 




slit structure 




hole structure 




SPP model 


^surf 




1.04 ±0.01 




1.05 ±0.01 




1.015 


Asurf (nm) 




819 ±8 




811 ±8 




844 


M 




0.13 ±0.01 




0.13 ±0.02 






K (/im) 




0.12 ±0.01 




0.020 ± 0.020 






V'intCTI-) 




0.81 ± 0.02 




0.55 ± 0.05 






TABLE II: Surface plasmon parameters 


Reference 




," 

tAg 




LabB(MHl) 


Lscat(Mni) 


Lrad(MHl) 


(7, 


-33.27 


1.31 


1.0154 


109 


2.56 X lO'* 


5.00 X 10" 


[14] 


-34. 


0.46 


1.015 


326 


2.6 X lO'' 


5.0 X 10" 


[15] 


-32.4 


1.74 


1.0158 


78.0 


2.51 X 10" 


2.68 X 10" 



"Measurements on silver films used in these experiments carried out at Caltech on a Sentech SE850 ellipsometer, 05 September 2005. 



with increasing distance out to about 3 — 4 /zm, but that the amplitude thereafter remains essentially constant. The 
solid curves in Figs. 6, 7 fit this amplitude decrease to an inverse distance dependence (Eqs. 8a, 8b). However, the 
damping might also plausibly fit an exponential decrease which would be expected from surface plasmon dissipative 
processes such as absorption by the silver film or scattering due to surface roughness. In order to check this possibility 
wc have measured the properties of the silver films used in these studies. Table II summarises these properties and 
compares them to previously reported measurements [14, 15]. The dielectric constant at 852 nm eAg = + ^Ag 
measured by ellipsometery and surface roughness parameters determined by atomic force microscopy (AFM). The 
root-mean-square (rms) height of the films was measured to be 5 = 1.29 nm and the correlation length a = 154.3 nm. 
From the imaginary term of the dielectric constant e^g and the parameters S, a the expected propagation lengths of 
surface plasmons against absorption, surface scattering, and reradiation, Lams, iscat, ^^rad, can be calculated [4]. It is 
clear from columns 5-7 of Table II that these loss processes cannot account for the observed damping within 3 /xm of 
hole-groove distance. 

The constant amplitude beyond ^3 — 4 /im is consistent with a persistant surface wave. Indeed we have recorded 
measurements (not presented here) of the surface wave persisting at least to ~ 30 /xm slit-groove distance. It is 
important to emphasise, however, that Asurf and risurf deviate significantly for those expected for a pure SPP on a 
plane silver surface. Intcrfcrometry measurements of the surface waves on "output side" slit-groove structures (not 
reported here) confirm the value of rigurf in Table I, and we believe that conventional, infinite-plane SPP theory [4] is 
not adequate to explain these results. We note that persistent surface waves over ~ 10 /um distances have also been 
reported in a double slit experiment [16] and interpreted as SPPs [17]. 

As indicated in Table I, intrinsic phase for the slit-groove and hole-groove structures respectively are (pf^j — 0.81 tt 
and <p(jj^ = 0.55 TT. Although one contribution to these phase shifts may be the CDEW "signature" phase shift of 7r/2, 
it is known from earlier studies that the specific form (width and depth) of the grooves themselves, can introduce 
phase shifts into the scattered wave [9]. We have determined the nature of these groove-induced phase shifts and 
resonances by measuring interference fringes arising from surface waves launched on the "output-side" of slit-groove 
structures. These results, that will be reported in a subsequent publication, support the existence of an intrinsic 
phase shift close to 7r/2. 

The interpretation that emerges from these results is that the subwavelength groove originates persistant, long- 
range surface waves by a two-step process: (1) the incoming TM polarised plane wave scatters from the groove and 
generates in its immediate vicinity on the surface a broad, CDEW-like distribution of diffracted evanescent waves, 
and (2) this broad-band local surface "emitter" excites, within a distance of ^ 3 — 4 /im, a long-range surface wave 
response. The near-term rapid amplitude decrease in the interference fringes of Figs. 6, 7 is evidence of this evanescent 
surface wave diffraction very near the groove. Persistant amplitude out to tens of microns is evidence for some kind 
of surface wave guided mode. It is significant to note that the wavelength and phase of the interference fringes do 
not shift over the entire range of the measurements. The initial diffracted surface wave components extend over a 
broad range of evanescent modes, kx > ko, including the conventional fcspp. Therefore it is to be expected that the 
local surface wave emitter excites this surface mode. We emphasise, however, that our measurements show that the 
wavelength of this persistant wave does not correspond to fcspp and that, when the phase lag associated the groove 
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itself is taken into account, the intrinsic phase of the surface wave with respect to the directly transmitted wave is 
close to 7r/2. The disaccord between Aspp and Xsurf is for the present a matter of speculation. Perhaps plasmon 
"leaky waves" [21] that transport energy very slowly away from the surface contribute to the spectrum of long-range 
surface excitation resulting in an effective wavelength shift; or perhaps, despite our ellipsometry measurements, the 
surface index of refraction of the metal film is slightly modified by some uncontrolled chemical or material process [22] . 
At a more practical level, these results indicate that it might be much easier to couple to surface guided waves than 
was previously thought. Conventional wisdom asserts that because the SPP lies to the right of the "light line" on 
the metal surface dispersion curve 5 ci j^i f iting or prism is needed to achieve efficient optical coupling [4]. A simple 
abrupt discontinuity in the surface, such as a slit or groove, appears to serve as an efficient coupler. Further studies 
are needed to understand the properties of the generated long-range persistant wave and to optimise the efficiency of 
this groove-coupling process. 

V. METHODS 
A. Structure fabrication 

The striictures consist of a single siibwavelength slit or hole flanked by one subwavelength groove. The grooves 
have a width of 100 nm and a nominal depth of 100 nm for the slit-groove structures and 70 nm for the hole-groove 
structures. The slit-groove distance {xgg) or hole-groove distance (xhg) is systematically incremented in the fabrication 
process. The subwavelength structures are fabricated by focused ion beam (FIB) milling (FEI Nova-600 Dual-Beam 
system, Ga+ ions, 30keV) into a layer of silver evaporated onto flat fused silica microscope slides. A low beam 
current (50 pA) was used in order to achieve surface features defined with a lateral precision on the order of 10 nm 
and characterised by near-vertical sidewalls and a minimal amount of edge rounding. Since it enables delivery of a 
variable ion dose to each pixel of the writing field, FIB milling conveniently allows the multiple-depth topography 
characteristic of the present devices to be formed in a single, self-aligned step. A 2-D matrix of structures is milled 
into the silver layer. Each matrix consists of 63 structures, nine columns, separated by 1.5 mm, and seven rows, 
separated by 2 mm. The first column contains only the slit with no flanking grooves. Light transmission through 
the slits in this column is used to normalise the transmission in the remaining columns. Variations in transmission 
through each of the elements in the "slits only" column provide a measure of the uniformity of the FIB fabrication 
process. Each entire matrix of structures is flanked on one side by a small round hole and on the other by a line 
grating for absolute reference positioning and angular alignment of the structure matrix with respect to the input 
laser beam. The square microscope slides themselves, commercially available from SPI Supplies, are 25 mm on a side 
and 1 mm thick. 



B. Measurement Setup 

Details of the experimental setup are as follows. Output from a diode laser source, temperature stabilised and 
frequency-locked to ^Si/2(F = 4) ^ ^P3/2(F = 4, 5) crossover feature in a Cs saturated absorption cell, is modulated 
at 850 Hz by a mechanical chopper, fed to a monomode optical fibre, focused and finally linearly polarised before 
impinging on the subwavelength structure mounted in the sample holder. The beam waist diameter and confocal 
parameter of the illuminating source are 300 fim and 33 cm, respectively. Throughout this series of measurements 
the laser power density was maintained ~ lWcm~^. The sample holder itself is fixed to a precision x-y translator, 
and multiple structures, FIB-milled in a 2-D array on a single substrate, are successively positioned at the laser beam 
waist. The optical response of the structures is synchronously detected by a photodiode and registered on a laboratory 
computer as indicated in Fig. 3. 
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